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^  Abstract 

The  liM«r  response  of  a  spatially  periodic  magnetized  Vlasov  plasma  distribution  function  is 
computed  to  second  order  in  the  electric  field.  The  results  for  a  specific  electric  field  are  then  compared 
with  the  results  of  computer  Simula  lion  for  different  amplitudes  of  the  electric  field.  Both  trapping  and 
resonant  heating  are  observed,  and  both  appear  to  contribute  (for  the  chosen  parameters)  to  limiting  the 
validity  of  linear  theory  at  larger  electric  field  amplitudes. 
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LINEAR  MAGNETIZED  PLASMA  RESPONSE 
TO  AN  OBLIQUE  ELECTROSTATIC  WAVE 


Introduction 

The  second  order  linear  response  of  an  nn magnetized  Vlasov  plasma  has  been  worked  out  [1].  It 
is  the  purpose  of  this  report  to  treat  the  magnetized  case.  This  extension  is  not  difficult  conceptually 
(provided  that  the  concepts  of  [1]  are  understood),  but  involves  quite  a  bit  more  algebra  than  is 
encountered  in  the  unmagnetized  case. 

Also  presented  here  is  simulation  work  showing  that  the  limits  of  linear  theory  are  imposed  by 
both  trapping  (as  in  the  unmagnetized  case),  and  by  perpendicular  heating.  Rough  estimates  are 
provided  for  the  wave  amplitudes  at  which  each  of  these  phenomena  become  important. 

Second  Order  Distribution  Function  for  Oblique  Electrostatic  Waves 

In  order  to  compare  the  details  of  Unear  theory  with  simulations,  the  second  order  perturbed 
distribution  function  /a(v, ,  vL )  is  needed  for  oblique  electrostatic  waves.  It  is  computed  in  essentially 
the  same  way  in  which  /3  was  calculated  for  unmagnetized  waves  in  [1].  As  before,  the  model  is 
1-d  periodic  with  no  collisions,  and  an  imposed  electric  field.  Now,  however,  there  is  a  magnetic 
field  at  an  angle  9  to  fc,  so  It  must  be  broken  up  into  a  kt  parallel  to  the  magnetic  field,  and  a  kL 
perpendicular  to  the  magnetic  field.  All  three  components  of  the  velocity  must  now  be  considered. 
For  the  purposes  of  analysis,  the  components  of  velocity  will  be  put  in  cylindrical  coordinates,  i.  e., 
v, ,  ,  and  0.  /2  contains  all  the  information  necessary  for  the  computation  of  linear  kinetic  energy 

and  mean  velocity.  The  second  order  perturbed  distribution  function  also  shows  how  the  kinetic 
energy  is  distributed  between  parallel  and  perpendicular  components. 

For  convenience,  the  coordinates  for  the  derivation  will  be  chosen  with  k  in  the  £  direction  with 
the  magnetic  field  at  angle  9  to  the  £  direction  in  the  x-z  plane.  The  phase  angle  0  will  be  defined 
so  that  when  the  magnetic  field  is  in  the  z  direction,  0  is  the  usual  angle  in  the  x-y  plane  with  0  =  0 
when  the  velocity  is  along  the  £  dirction,  and  0  =  ir/2  when  the  velocity  is  in  the  y  direction  (0  is 
90*  out  of  phase  with  the  phase  with  respect  to  position). 
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The  result  is  obtained  by  the  method  of  characteristics,  just  as  in  the  case  of  unmagnetized 
waves,  although  now  the  characteristics  are  helices  rather  than  straight  lines.  The  final  result  is 
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where  the  bar  denotes  averaging  over  both  position  and  phase.  (See  Appendix  A  for  the  complete 


calculation.) 


Note  that  as  in  the  case  of  the  unmagnetized  wave,  the  operator  on  /o  is  like  a  diffusion  operator. 
In  fact,  one  could  make  a  quasilinear  diffusion  equation  from  the  equation  for  dfa/dt  by  substituting 
f  lot  ft  and  /o-  The  form  of  the  operator  implies  that  all  else  being  equal,  the  diffusion  in  the  vx 


direction  is  greater  for  larger  n. 


Since  Ji  is  now  a  function  of  two  variables  (v,  and  vx),  it  will  not  be  so  useful  a  diagnostic 
in  particle  simulations  due  to  the  noise  inherent  in  trying  to  fill  a  two-dimensional  space  with 
only  130,000  particles.  Scalar  quantities  such  as  the  mean  parallel  velocity,  end  the  parallel  and 


perpendicular  kinetic  energies  (the  mean  perpendicular  velocity  is,  of  course  zero  in  a  magnetized 
plasma),  are  much  better  suited  to  be  used  as  diagnostics.  Fortunately,  they  are  easy  to  compute 
from  the  distribution  function.  Assuming  a  single-temperature  Maxwellian  distribution  for  /o,  and 
letting  ft  =  kj_vt/ue,  u  =  v%/vt,  z  —  u/knvt>  w  s  ue/ktvt  and  a  =  y/ktvt  the  results  are 
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These  integrals  are  well-behaved,  and  can  be  evaluated  numerically  with  relative  ease. 


As  with  unmagnetized  waves,  a  formula  can  be  derived  for  /2(f  —  oo).  The  derivation  follows 
that  of  /2.  The  result  is 
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where  E  is  the  Fourier  transform  of  the  electric  field  envelope  E. 

This  last  expression  is  much  better  suited  to  numerics!  comparison,  as  the  particle  distribution 
function  is  difficult  to  obtain  accurately  at  any  given  instant.  As  with  the  instantaneous  formulae  for 
and  £x,  the  integrals  for  these  quantities  as  t  — ►  oo  are  relatively  easy  to  evaluate  numerically 
for  the  electric  field  envelope  which  will  be  used  in  the  simulations. 

It  is  again  worth  noting  the  resemblance  of  this  result  to  those  of  quasi- linear  theory  (see,  for 
instance,  Kennel  and  Engelmann  [2]). 

Theory  for  Trapping  and  Perpendicular  Heating  Due  to 
Oblique  Electrostatic  Waves 

Trapping  can  occur  in  waves  in  a  magnetized  plasma,  such  as  oblique  electrostatic  waves,  as  well 
as  in  unmagnetized  waves.  Another  phenomenon,  perpendicular  heating,  is  unique  to  the  magnetized 
case,  and  must  also  be  considered.  It  is  the  primary  goal  of  this  research  to  find  the  field  strength  at 
which  trapping  and  perpendicular  heating  become  important.  First  trapping  will  be  analyzed,  then 
perpendicular  heating. 

The  basic  concept  behind  trapping  is  the  same  in  the  magnetized  case  as  in  the  unmagnetized. 
The  mechanism  for  a  magnetized  plasma  is,  however,  somewhat  more  complicated.  Instead  of  a 
single  resonant  velocity,  there  are  an  infinity  of  them,  one  corresponding  to  each  harmonic  of  the 
cyclotron  frequency.  All  these  resonant  velocities  satisfy  w  -  kt  v„  -  nuc  =  0.  In  the  rest  frames  of 
particles  at  each  of  these  resonant  velocities  the  wave  has  the  same  apparent  k,  but  an  apparent 
frequency  of  nwe  where  n  is  the  harmonic  number  (in  the  non-relativistic  approximation,  a  boost 
parallel  to  the  magnetic  field  does  not  alter  the  fields).  To  understand  the  trapping  in  the  parallel 
direction,  it  is  necessary  to  average  the  force  on  a  trapped  particle  over  the  short  time  scale  (the 
gyromotion).  From 

£v-s1(e  +  v-xb)  (6) 


the  equation  in  the  parallel  direction  is  obtained  by  dotting  with  the  unit  vector  in  the  direction  of 
the  magnetic  field,  6: 


Prom  this  point  on,  the  coordinates  will  be  chosen  so  that  £  is  in  the  z  direction,  and  the  electric 


field  is  in  the  x-z  plane.  Thus,  k-z  **  j fe,z(t)  +  fcxx(t).  Since  we  are  at  the  resonant  velocity  looking 
at  particles  which  are  at  or  near  resonance,  z(t)  (the  z  coordinate  of  the  particle  in  the  resonance 
frame)  is  slowly  varying  (Le.,  v,  is  small)  since  it  is  parallel  to  the  magnetic  field,  and  only  x(t)  and 
y(t)  will  be  rapidly  varying.  Only  x(t)  is  of  interest,  and  setting  x(t)  =  A'(t)  +  ^  sm(u/et  4-  4>)  where 
X(t)  is  the  guiding  center  position  and  therefore  slowly  varying,  gives 

E(t)  =» 
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Averaging  this  over  a  gyroperiod  (and  taking  the  real  part)  gives 


E  -  Ecmiktz  +  kxX  +  rut>)Jn 


Thus, 


£«,  -  ±EaM9a»(ktZ  +  KX  +  n*)Jn 


and  for  the  most  deeply  trapped  particles, 


where  z'  is  the  z -component  relative  to  the  bottom  of  the  trapping  well.  This  is  the  same  formula 


as  for  one  dimension,  except  for  the  factors  of  coed  and  Jn.  These  factors  are  easy  to  understand 


in  physical  terms.  The  cm  8  factor  exists  because  of  the  angle  between  the  electric  field  and  the 


direction  of  motion.  It  is  the  component  of  the  electric  field  along  the  parallel  direction  instead  of 
the  entire  electric  field.  The  second  factor  results  from  the  averaging  of  the  electric  field  over  the 


particle  orbit,  which  may  be  comparable  to  a  wavelength  in  radius. 


The  same  quantities  Av(r  and  A4(r  (the  maximum  trapping  width  and  the  cumulative  trapping 
phase  shift)  can  be  extended  directly  from  the  one  dimensional  case  to  this  equation  for  v, .  By  direct 


analogy. 
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As  in  the  one  dimensional  case,  if  either  At Hr  or  A$tr  is  small  enough,  the  distribution  will  teem  to 
follow  the  linear  theory  closely,  but  in  the  self-consistent  case,  A$(P  should  be  the  only  important 
parameter. 

As  mentioned  in  the  beginning  of  this  section,  the  magnetic  field  introduces  a  new  effect: 
resonant  perpendicular  heating  (or  cooling).  The  derivation  is  similar  to  the  derivation  for  the 
trapping,  beginning  with  the  same  equation  of  motion  (6),  but  dotting  it  with  vx  instead  of  6.  The 
result  is 


=  —E sm0vJe<(*«*+*J-x+n*)  cos(wet  +  (14) 

xn  ' 

Averaging  again  over  a  gyro  period  gives 

5  (ft)  =  ^EemOcosiktZ  +  k^X  +  Tut,)^,,  (15) 

This  formula  has  some  interesting  ramifications.  For  particles  which  are  not  trapped,  z  will 
vary  fairly  rapidly,  and  the  average  of  the  right-hand  side  will  be  close  to  zero.  Particles  which 
are  trapped,  however,  have  a  limit  to  the  z  coordinate,  and  so  the  right-hand  side  may  have  a 
substantial  average.  This  formula  is  also  interesting  in  the  extreme  non-linear  case  because  of  the  Jn 
factor.  Since  the  argument  of  the  cosine  is  nearly  constant  for  given  particles,  it  will  either  gain  or 
lose  energy  monotonically  until  the  value  of  vx  approaches  a  value  such  that  Jn(Jfciux/u;<;)  vanishes 
(except,  of  course,  for  the  case  n  =  0).  Since  the  trapping  in  v„  also  has  a  factor  of  J„  in  it  (Equation 
(11)),  the  effect  of  the  resonance  in  the  perpendicular  direction  is  to  reduce  the  effects  of  trapping 
in  the  parallel  direction,  except  for  the  important  case  of  n  =  0.  Since  the  effect  of  this  equation  is 
to  reduce  the  effects  of  trapping,  yet  it  relies  on  trapping  for  its  effect,  it  is  somewhat  self-limiting. 
The  following  derivation  of  a  condition  for  linearity  does  not  take  this  into  account,  and  so  should 
be  quite  conservative.  At  first  inspection,  it  might  also  seem  that  the  perpendicular  damping  would 
not  disappear  when  the  angle  of  propagation  is  perpendicular  to  the  direction  of  the  magnetic  field. 
It  is  assumed,  however,  in  deriving  (15)  that  the  particle  being  followed  is  in  resonance  with  the 
wave,  and  since  is  zero  for  perpendicular  propagation,  resonance  requires  that  w  -  nuie  =  0.  This 
is  never  true  for  Bernstein  modes  which  are  propagating  perpendicular  to  the  magnetic  field,  so  the 
mode  will  remain  undamped  as  theory  predicts. 


The  perpendicular  resonance  ran  be  said  to  have  become  non-linear  when  the  perpendicular 
velocity  has  changed  by  enough  to  significantly  alter  the  value  of  Jn(kxv, or  if  the  particle 
picks  up  perpendicular  energy  comparable  to  the  thermal  energy,  i.e.  the  distribution  function  is 
significantly  altered  on  a  scale  large  enough  to  affect  the  bulk  plasma,  and  thereby  the  dispersion 
relation.  For  a  rule  of  thumb,  the  factor  of  cos^r  +  kxX  +  n4>)  can  be  o  mi  ted,  and  the  factor  of 
Jn  can  be  replaced  by  1/2  (since  the  Jo  term  does  not  contribute),  leaving 

(16) 

dt  TO 

The  solution,  plugging  in  E  **  Eq  exp(7<)  is, 
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For  linear  theory  to  hold,  this  quantity  should  be  much  less  than  1,  but  it  may  hold  reasonably  well 
for  larger  values  since  all  that  is  really  required  is  that  Avx  <C  which  is  a  weaker  condition 

than  (18).  It  should  also  be  noted  that  the  higher  values  of  n  are  not  as  important,  so  for  the  rule 
of  thumb,  n  can  be  omitted  from  this  formula. 

The  second  condition,  that  the  change  in  energy  be  much  less  than  the  thermal  energy,  is  derived 


trivially  from  the  same  formula: 


(vx\2  qweEosini 


Again,  n  may  be  omitted  for  a  rule  of  thumb.  The  final  result  is  that  if 


qEo*c  _.2\ 
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and  the  trapping  condition  A4<P  ir  is  satisfied,  then  linear  theory  will  be  accurate. 

Note  that  these  formulae  depend  on  £0/7  rather  than  on  i/Eo/y  as  with  the  parallel  trapping. 
Thus  for  «rnall  damping  rates  and  wnall  field  amplitudes,  the  perpendicular  heating  will  cause  non¬ 
linear  effects  before  trapping  can  set  in. 


Numerical  Particle  Simulations  of  Oblique  Electrostatic  Wave 


The  sort  of  simulation  as  was  done  for  unmagnetized  waves  in  [1]  can  also  be  done  for 
oblique  electrostatic  waves.  As  with  the  unmagnetized  waves,  a  large  number  of  particles  is  necessary; 
but,  unlike  the  unmagnetized  simulation,  the  distribution  function  is  in  two  variables  -  v,  and  vx. 
As  before,  the  field  amplitude  at  which  non-linear  effects  set  in  is  of  interest.  The  situation  is  not 
so  simple  as  for  the  unmagnetized  wave,  and  so  simulation  should  be  a  useful  adjunct  to  theory. 

Again,  the  simulation  model  follows  the  theoretical  model  (as  described  in  the  previous  section. 
It  is  1-d  (with  all  three  velocity  components,  however),  collisionless,  and  periodic,  with  a  neutral¬ 
izing  background  charge  density  representing  immobile  ions.  The  initial  distribution  function  /o  is 
Maxwellian,  and  the  imposed  electric  field  is  again  E  —  Eo  exp(— 7|t|)  exp (i(kx  —  wt)). 

The  distribution  function  cannot  be  as  easily  plotted  and  interpreted  in  this  case  as  it  was  in  the 
case  of  unmagnetized  waves,  because  it  is  a  function  of  both  t/(  and  vx .  This  problem  is  compounded 
by  the  fact  that  the  number  of  particles  is  not  enough  to  smoothly  fill  the  two-dimensional  velocity 
space.  (Two  examples  of  contour  plots  of  f(vt,  vx)  will  be  shown,  but,  as  will  be  apparent,  they  are 
not  as  informative  as  the  distribution  function  plots  were  for  the  unmagnetized  wave.)  The  mean 
parallel  velocity  (which,  of  course,  represents  direct  current  drive),  the  parallel  kinetic  energy,  and 
the  perpendicular  kinetic  energy  (which  <*an  generate  current  indirectly  when  collisions  are  present) 
are,  however,  easily  calculated  and  just  as  informative  as  were  the  mean  velocity  and  kinetic  energy 
for  the  unmagnetized  wave.  Each  of  these  can  be  calculated  from  the  theory  (through  through  the 
numerical  evaluation  of  some  well- behaved  integrals)  and  horn  the  simulation  both  as  the  imposed 
w  ave  is  growing  exponentially,  and  as  t  —  oo. 

The  parameters  chosen  for  these  runs  are:  t/u/e  =  1,  &xv*/u/e  =  2.4  (so  kvt/u/e  =  2.6  and 

0  as  67°),  ui/u)c  =  1.5,  ~f/uic  =  -25.  Five  runs  were  made  at  differing  values  of  Eo-  These  values 
were  chosen  such  that  qEo/mvttje  would  be  equal  to  0.125,  0.25,  0.5,  1,  and  2  (note  that  this  is 
a  slightly  different  normalization  from  that  of  the  unmagnetized  wave  simulations  —  but  that  the 
two  can  be  compared  through  the  formula  for  A$(r).  As  with  the  unmagnetized  simulations,  the 
parameters  were  chosen  for  ease  of  simulation  rather  than  for  realism,  and  bear  no  resemblance  to 
any  self-consistent  wave. 

A  least  squares  fit  to  a  growing  exponential  with  exponent  7  was  made  on  the  last  e-folding  in 
E  of  the  growing  part  of  each  the  three  diagnostics.  The  purpose  of  these  curves  was  to  obtain  a 
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more  accurate  value  for  the  diagnostics  (which  they  did)  and  to  provide  an  estimate  of  the  er^r  due 
to  noise  which  could  be  used  as  an  estimate  of  the  error  of  the  diagnostics  as  t  — *  oo  (which  they  did 
not).  The  standard  deviations  of  the  diagnostics  about  the  fit  exponential  curve  should  be  a  fairly 
good  measure  of  the  error  of  any  given  point  on  the  diagnostic  curves.  The  instantaneous  value  of 
the  diagnostics  is  actually  plotted  in  Figure  1,  rather  than  the  value  obtained  from  the  least  squares 
fit,  for  two  reasons:  first,  the  assumption  under  which  it  was  deemed  appropriate  in  the  first  place 

—  that  a  fixed  level  of  noise  was  superimposed  on  a  basically  correct  curve  —  appears  to  be  false, 
and  second,  the  correct  curve  is  expected  to  deviate  from  a  growing  exponential,  and  thus  there 
would  be  a  systematic  error  in  the  least  square  value. 

Figure  1  shows  the  values  of  the  diagnostics  obtained  from  simulation  (appropriately  normalized) 
and  the  theoretical  linear  values.  All  the  diagnostics  converge  toward  the  linear  values  quite  well 
as  the  field  amplitude  decreases,  except  for  the  parallel  energy,  which  seems  to  show  either  an 
oscillation  with  field  amplitude,  or  fairly  strong  noise.  Since  another  run  with  a  slightly  modified 
particle  loader  (bit-reversed  in  the  x-direction  instead  of  the  direction  parallel  to  the  magnetic  field 

—  see  the  technical  note  at  the  end  of  this  section)  yielded  values  which  deviated  with  similar 
amplitude  but  without  sign  of  oscillation,  it  seems  likely  that  the  error  is  due  to  noise  caused  by 
some  small  systematic  error  in  the  loading  scheme.  This  is  initially  somewhat  surprising,  given  the 
large  number  of  particles  in  the  system,  but  quite  reasonable  when  it  is  remembered  that  only  a 
small  fraction  of  those  particles  are  resonant  with  the  wave.  Note  also  that  the  perpendicular  kinetic 
energy  decreases  with  increasing  field  amplitude.  Since  perpendicular  heating  can  drive  a  current 
once  the  effects  of  collisions  are  taken  into  account,  this  loss  of  efficiency  is  important.  More  will  be 
said  about  the  results  shown  in  Figure  1. 

Figure  2  shows  the  parallel  kinetic  energy  as  a  function  of  time,  along  with  the  exponential  curve 
which  was  fit  to  the  parallel  kinetic  energy  using  least  squares,  for  the  lowest  amplitude  excitation 
over  the  last  e-folding  of  the  electric  field  before  t  =  0.  The  estimated  error  from  this  least  squares 
fit  is  ~3%,  which  is  in  fair  agreement  with  the  deviation  of  this  diagnostic  from  theory  at  (  =  0 
(—6%),  but  far  smaller  than  the  deviation  from  theory  at  t  =  oo.  That  some  low  frequency  noise  is 
present  is  plainly  visible.  That  it  is  only  low  frequency  and  that  it  exists  at  all  in  what  should  be  a 
uniform  system  is  more  evidence  for  the  noise  being  the  result  of  systematic  error. 

Figure  3  shows  contour  plots  of  f(v^,vx)  before  the  simulation  and  at  t  —  oo  for  the  highest 
excitation  amplitude,  dearly  showing  the  effects  of  the  wave.  Some  current  drive  effect  is  also  clearly 
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visible. 

When  Figure  1  is  compared  with  the  results  of  the  unmagnetized  simulation,  it  is  seen  that, 
as  with  the  unmagnetized  simulation,  the  linear  theory  holds  up  to  a  certain  field  amplitude,  then 
the  results  start  diverging.  Applying  Formula  (13)  to  the  simulation,  it  is  found  (setting  Jn  ~  .5) 
that  A4|,  ranges  from  1.67  to  6.67,  and  that  comparable  deviation  is  found  at  comparable  values 
of  This  is  a  strong  confirmation  of  the  validity  of  Formula  (13). 

The  perpendicular  heating  condition  (20)  can  aslo  be  evaluated.  The  values  of  the  left-hand 
side  range  from  0.19  to  3.08,  and  the  values  of  the  right-hand  side  (which  should  both  be  much 
larger  than  the  value  of  the  left-hand  side  for  linearity  to  be  guaranteed)  are  0.17  and  1.0.  The 
good  agreement  of  the  simulation  with  linear  theory  indicates  that  either  an  error  was  made  in  the 
derivation  of  the  condition  (20),  or  that  the  many  approximations  made  in  deriving  (20)  were  very 
conservative. 

A  mildly  surprising  feature  of  the  simulation  is  that  the  parallel  energy  increases  with  increasing 
non-linearity.  This  is  easily  explained  when  it  is  realized  that  the  distribution  function  is  actually 
made  steeper  at  some  resonances  according  to  linear  theory.  When  trapping  sets  in,  these  places  do 
not  steepen  as  linear  theory  predicts,  and  may  even  flatten  out  instead  of  steepening.  Thus  these 
resonances  contribute  more  heating  than  expected.  This  effect  does  not  contribute  to  the  current 
drive,  as  the  resonances  which  ought  to  steepen  the  distribution  function  are  primarily  at  velocities 
opposed  to  the  direction  of  current  drive. 

Technical  note: 

Because  of  the  high  noise  level  associated  with  random  loading  in  this  number  of  phase  space 
dimensions,  the  particles  were  loaded  with  some  care.  The  velocities  in  the  parallel  direction  (rather 
than  the  x  direction)  were  loaded  in  bit-reversed  order,  the  phase  of  the  perpendicular  velocity  was 
chosen  in  3- re  versed  order,  and  the  magnitude  of  the  perpendicular  velocity  was  chosen  in  5-reversed 
order.  This  arrangement  should  introduce  the  minimum  noise  in  the  parallel  direction.  (The  bit- 
reversed  loader  has  fewer  problems  with  recurrence  than  the  3-  and  5-  reversed  loaders.)  See  Birdsall 


and  Langdon  [3]  for  details  on  bit-reversing. 
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Summary  and  Conclusions 


Linear  theory  has  been  used  to  compute  the  second  order  perturbed  distribution  function  for 
an  oblique  electrostatic  electron  wave  both  for  an  exponentially  growing  wave  (Equation  (1))  and 
for  a  wave  pulse  of  arbitrary  shape  (Formula  (5)).  The  results  are:  for  the  growing  wave 
t  /  v  v-''  ( .  d  nue  d  \ 

/aK.vJ  -  -y  •  +  ‘ 

1 - “c (^TiT  +  n“~  fH  )  fo(vi,Vj.)  (21) 

(w,  ~nue  -Je,i;l)2  +  7a  \  "  dvt  vx  Qv^J 

and  for  the  wave  pulse 

1  _  1  f  V  (fr  9  9  \  . 

3  4  kPm2  "  \  1  du,  «j.  / 

[j2  (^)  |£(«  -  -  *.».)!’  (*,£  +  /.(«•.».)]  («) 

This  linear  theory  has  been  tested  via  particle  simulation  on  its  predictions  regarding  several 
second  order  quantities.  The  electric  field  varied  in  time,  though  it  was  uniform  in  space.  The 
results  were  in  agreement  with  theory,  although  the  tests  have  not  been  self-consistent  in  the  sense 
that  the  electric  field  was  imposed  rather  than  solved  for  using  Poisson’s  equation.  The  onset  of 
non-linearity  was  observed,  and  the  threshold  level  fit  the  predictions  of  a  simple  trapping  model. 
This  model  yields  two  parameters  which  predict  the  behavior  of  the  plasma  in  the  presence  of  oblique 


electrostatic  waves:  the  cumulative  trapping  phase  shift 

=  -\  —k^EocoadJn 

yy  m 

and  the  perpendicular  energy  change 


2  1  qE0  TIUr. 

A(<) - 2 —  •  -r- 

7  m  Ac 


7  |  m  |  Ac 

If  both  of  these  parameters  are  small,  the  behavior  should  be  linear,  and  indeed  the  results  are 
consistent  with  this  prediction. 

It  is  fair  to  say  that  the  one  dimensional  effects  causing  non-linearity  are  now  understood.  Un¬ 
fortunately,  these  conclusions  cannot  be  blithely  applied  to  three  dimensions.  While  the  cumulative 
trapping  phase  shift  and  perpendicular  energy  change  are  functions  only  of  the  electric  field  strength 


•▼V 


11 


Appendix  A:  Second  Order  Perturbed  Distribution  Function  for  Oblique 
Electrostatic  Wave 


For  convenience,  the  coordinates  for  the  derivation  will  be  chosen  with  k  in  the  z  direction  with 
the  magnetic  field  at  angle  9  to  the  z  direction  in  the  x-z  plane.  The  phase  angle  0  will  be  defined 
so  that  when  the  magnetic  field  is  in  the  z  direction,  0  is  the  usual  angle  in  the  z-y  plane  with 
0  =  0  when  the  velocity  is  along  the  z  dirction,  and  0  —  r/2  when  the  velocity  is  in  the  y  direction 
(0  is  90*  out  of  phase  with  the  phase  with  respect  to  position).  The  unperturbed  orbit  in  such  a 
coordinate  system  is  such  that 

t>,(0  as  vxcos(u >e(tf  —t)  +  0(t))sin0  +  v,  cos#  (ol) 

z(0  =  z(t)  +  [sin(ure(t'  -t)  +  0(t))  -  sin0(t)]  sin 9  +  -  t)cos0  (a2) 

0(0  «0(t)  +we(t'-t)  (o3) 

with  v,  and  vx  constants  (recall  that  He,  =>  kccaO  and  s  ksinO). 


The  first  equation  to  be  solved  is 


4ft  _  _2_b  £/o 

dt  m  9v 


E  =  Ekei(h—Wt) 


In  cylindrical  velocity  coordinates,  the  differential  operator  becomes 

-  d  1  C,  a  3  L  .  Id]  , 

+  <a6) 

Plugging  this  operator  into  the  equation  for  f\  and  assuming  that  /0  does  not  depend  on  0,  yields 

T*  =  ““  T  [*•  ST  +  cos0*x  /o  •  («7) 

dt  m  k  [  dvt  9vx 

Plugging  in  x(f)  and  0(0  and  solving  by  the  method  of  characteristics, 

m  k  \  1  dv,  J.m 

+  k^  J‘  cos(0  +  we(t' 


__!L 

m  Jk 


x 
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/t  df*  e*(fcJ.*x  rn(*-«.r)-k,.,f+«r]  fa. 

X  &>%  Jo 

+  fcx  J~  cas{4>  -  ucr)et^hj-9-L  n'+^l  <frj 


_±  “  x 
m  fc 


E  [*■§£■ +  ^  IS]  *  (^r) -f  — ‘-,T'fr 

-  [*£ + tt^t]  *  (a8, 

Note  that  it  «u  necessary  to  assume  Imu  >  0,  i.e.  a  growing  wave.  It  will  also  be  useful  when 
calculating  f2  to  write  f\  as 

/, .  -i-i  £«•<*•-'*> 


EE 


r  a  no.,  a 


u>  —  nu/e  —  At|»| 


L*  +«:  «1A 

[  dvt  I/x  OVj. 


Note  that  the  first  order  perturbed  distribution  function  is  a  function  of  the  velocity  phase  angle, 
which  the  zero  order  distribution  function  is  not. 

The  equation  for  the  second  order  perturbed  distribution  function  is 

dt  m  &vn 

*  -^£Ree<(*»'“'‘)  •  Re  Jx  (alO) 

The  real  part  of  all  quantities  being  previously  implied. 

Before  integrating  the  equation,  it  is  best  to  average  over  x ,  as  the  spatial  structure  of  /3  is  of 
no  interest.  It  is  also  assumed  at  this  point  that  k  is  real.  Now 

m  2k  I  V  dv )  ^  ui  -  nue  -  fc,  r, 

\  fl  fft 

where  7  =  Imw.  The  velocity  phase  angle  0  can  also  be  integrated  over  now  since  the  phase 
structure  of  f2  is  of  no  value  in  computing  the  quantities  of  interest.  (One  must  actually  be  careful 


to  denote  4>(t)  as  4(to)  +  u ;e(t  —  to)  and  average  over  ^(to),  bat  the  end  resalt  is  the  same.)  The 

differential  operator  can  be  rewritten  as 

?  d  _  L  d  ,  ,  x  d  fcj.sin<£  3 

*  •  3—  +  fcj.  C06<£— - — 

&V  (TV,  3vx  Vj.  3<p 

.  5  ,  fcj.  3  .  hi  a  .  .  , 

(ol2) 

Integrating  over  4>  will  eliminate  the  d/Q4>  term  so  that  only  two  integrals  mast  be  calculated: 

e*(-”)*g  =  6n.m  (al3) 

and 

cos^e*( - >♦  g  =  +  Wi)  (°14) 

Using  these,  the  sums  over  m  can  be  evaluated: 

£  Jm<5n.m  =  in  (al5) 

m 

Er  if  \  _  *^n+l(2)  ~f~ 

2  1  +  ®n,m+ll - 2 


Plugging  these  in, 


d/a  .£?  g2  W, 
dt  2  IAnJ  \  1 3v,  Vj.  3vj.  / 

(  A  )  f  l  &  Wc  9  \  2-r< 

- t —  (  5T  + - AT"  /o  x 

u t  -  n«c  -  fc,v,  \  3v,  Vj.  av j.  / 


Now  taking  the  real  part, 

*4 Wt  A^u^A^i 

dt  4  fc2m2  \  1  3v „  Vj.  3vj.  / 


- 1 

aS 

1 

*  A. 

^  3  \ 

|^(w,  -  nwe  -  A,v,)2  +  72  \ 

vj.  3vxJ 

and  so 


x  .  .  E*  <j*  /,  3  nu/e  3  N 

/2(wa,vJ  =  + 

(*t^)  3  _mje  d  \  tl 

(uir  -nu/e-  fc,v,)2  +  72  V  "  3v„  Vj.  3vj.  /  '°  U"’ 1,1 


The  kinetic  energy  of  the  wave  can  now  be  calculated  in  order  to  compare  the  result  with  that 
of  the  energy  derived  from  the  linear  dispersion  relation.  The  kinetic  energy  can  be  separated  into 
parallel  and  perpendicular  components.  Starting  with  the  parallel  component, 


£,  =  JJ  ^mv2f2  dv,  vxdvj_ 


Breaking  up  ]2  into  two  terms,  the  d/&v±  term  can  be  integrated  immediately,  yielding  zero  since 
n*A>(0)  =  0  for  all  n.  Integrating  the  d/dvt  term  by  parts  and  substituting 


( ,  Q  nu>c  3  \  1 

[kt9^  +  fo  "  +  «*)/o 


yields 


Setting 


5"  "  4  *  A*  \>]J  +  \~)  f*dV' 


V^ir8  vf 


separates  the  integrals: 


Setting  m  =  k^vt/uie  and  y  =  reduces  the  first  integral  to 

/ ■*"  (^r)  ^  ~  /  Jn(i*y)e~lr  ydy 

=  (a25) 

The  algebra  for  the  second  integral  can  be  simplified  if  some  new  symbols  are  introduced.  Let 

u  =  vjvt,  z  =  w/fc„v,,  u>  =  ve/ktvt  and  a  =  7/fc, ut  then  the  integral  becomes 
1  f  u(t*  4-  nu>) 
y/2ir  J  («  +  nui  -  z)2  4-  at2* 

1  f  (2z  —  nt»)(tt  —  r  +  nu)  4-  z(z  -  ntt/)  —  ar2  _«i 

=  1+7si - [.-;+«.)»+«■ - *  * 

=  1  4-  (nt»  -  2z)  Re  — -L=  f - - - —  du 

V27r  J  u-:+niD-ia 

+  1  r - 

a  V  2?r  y  u  -  z  4-  nw  —  »a 

=  1  4-  (nw  -  2 z)  Re  lV(z  -  nw  4-  ia)  +  — — — —  lmN(z  -  nw  +  ia)  (o26) 

a 


IS 


Taking  the  limit  as  a  — *  0  for  comparison  later  with  the  results  of  the  calculation  using  the  linear 
dielectric  function, 

Im  JV(z  —  nw  +  tor)  — ►  Im  [lV(z  —  nw)  +  iaN'(z  —  nwj] 

— *  Im  [lV(z  —  nw)  —  ia(l  +  (z  —  nw)N(z  -  nw))] 

— *  ImiV(z  —  nw)  -  a(l  +  (z  -  nw)  Re N(z  -  nw))  (“27) 


(where  N(z)  =  l/y/2Z(z/V2)  and  Z  is  the  standard  Plasma  Dispersion  Function)  and 

£’•  —  -*(«-««)  +  [2*  -nw)3]  N{z  -nw)+ 

- >.1  (a28) 

The  perpendicular  energy  is  evaluated  in  mnch  the  same  way: 

A.  =  jj  do,  VidOj.  d* 

rf  nuieJ”  (^)  A  d  wt  fl\  ^  . 

4  k2  ^JJ  (w  -  nw«  -  fc,w,)3  + -y3  (  +  vx  dvj.)  *  0  ,Vx  x 

~  4  kaX2o  >  J  {u-z  +  nw)*+a*V2;e 


E2  1 
4  A^A3 

£?  1 
4  fc3  A3 


y  In(fit)e~l,i  |nwReiV(z  -  nw  +  ia)  +  ImJV(z  —  nw  +  »a)  j  (“29) 
*1 

£in(MJ)e-** 


x  |-znw  +  [nw  -  znw(z  -  nw)]  Re.V(z  —  nw)  +  \j^~~ -e_^(,-n,•,,  | 


(“30) 


The  total  energy  is  thus 
£  =  4*  fi 


A,  a  ,  nwe  a  \  , 

4  (w-  nwe  -*lt',)3  +  r3  \  »1  0Vj./  0  ,Vi  Vi 

1  ^  /7(*,v,+nwe)»./*(i 

4  ^A|?yy  (w-nwe-fc,t',)3  +  73  *  °  V|  Vi  A 


|l  +  2zRbN(z  -  nw  +  ia)  +  *  •  — —  ImlV(z  -  nw  +  iq)|  (<*31) 

|l  -  z3  +  [2z  -  z*(z  -  nw)]ReJV(z  -  nw)  +  (<*32) 

To  ensure  that  the  calculation  of  /j  is  correct,  the  kinetic  energy  of  a  electrostatic  wave  will 
now  be  calculated  using  the  dielectric  function: 

^  fJ-jO*/  Ree)  -  «oj  («*33) 


Wi* mdmu  =  Ul  Im  e  ■ 


r*  £(oa 


where  is  the  wave  kinetic  energy  density  and  Wc,n^T«  is  the  energy  which  has  been  absorbed 


by  the  distribution  function  due  to  Landau  damping  The 


that  ail  changes  in 


the  electric  field  amplitude  are  infinitesimally  slow,  and  the  result  is  only  the  total  energy,  rather 
than  h  itself.  The  result  will  be  a  completely  independent  check,  however,  on  the  energy  as  derived 
from  the  second  order  perturbed  distribution  function.  The  result  is 

-?sAr{'- m -W  '(t?) 


where 


=  V2Z(ti) 


and  Z  is  the  standard  Plasma  Dispersion  Function.  The  contour  of  integration  goes  under  the  pole 
at  C  =  z  as  with  the  Z  function.  Use  has  also  been  made  of  the  identities 

N'(z)  =  -(1  +  zN(z))  (o37) 


5] /«(*)€-  -1 


To  amplify  the  notation,  let  z  *  w/fc, vt,  w  a  ue/ktvt,  and  fi  a  k±vt/u>c,  then 

& 


W, 


T*51ij  j1  “  *3  +  D2*  ~  **<*  "  nw)W{z  ~  |  («39) 


0/ 

«  j  *  -Im« 


The  purely  kinetic  (Landau  damping)  part  of  the  energ'  for  an  exponentially  growing  wave  is 

(o40) 

v  7 

where  7  is  the  growth  rate.  For  real  argument, 

InxiV(z)  -  Jle-*  (o41) 


so,  setting  a  a  7/fc ,u«, 

"*•— -t  T  v  ''-Urr  • 

**  1 


4  a 


(o42) 


Equations  (a39)  and  (a42)  sum  to  Equation  (a32),  thus  verifying  it. 
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